Hot and dense gas of quark quasi-particles 
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Some features of hot and dense gas of quarks which are considered as the quasi-particles of the 
model Hamiltonian with four-fermion interaction are studied. Being adapted to the Nambu-Jona- 
Lasinio model this approach allows us to accommodate a phase transition similar to the nuclear 
liquid-gas one at the proper scale. It allows us to argue an existence of the mixed phase of vacuum 
and baryonic matter (even at zero temperature) as a plausible scenario of chiral symmetry (partial) 
restoration. 
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Understanding in full and describing dependably the 
critical phenomena (chiral and deconfinement phase tran- 
sitions) in QCD is still elusive because of a necessity to 
have the corresponding efficient non-perturbative meth- 
ods for strongly coupled regime to be analyzed. For 
the time being such studies are pursued by invoking di- 
verse effective models. The Nambu-Jona-Lazinio(NJL)- 
type models are certainly playing the most advanced role 
in this analysis [1]. This approach deals with the four- 
fermion interactions in lieu of a gluon field QCD dynam- 
ics and does not incorporate (albeit being supplemented 
by the Polyakov loops it does [2|) the property of con- 
finement. In the meantime it is quite successful in realiz- 
ing the spontaneous breakdown of chiral symmetry and 
its restoration at non-zero temperatures or quark densi- 
ties. Apparently, the non-renormalizable nature of the 
NJL model supposes another approximations introduced 
to be solved and might lead to some conclusions which 
are dependent on the regularization scenario sometimes. 
Hence, it requires a steadfast control of all inputs done 
to have the consistent physical results and to avoid the 
reasons for skeptical attitude. 

Instructive example was given in Refs. |1 and Q m 
which the ground state of model Hamiltonian with four- 
fermion interactions was studied in detail. The quarks 
were treated as the quasi-particles of this Hamiltonian 
and unexpected singularity (discontinuity) of the mean 
energy functional as a function of the current quark mass 
was found out. In particular case of the NJL model new 
solution branches of the equation for dynamical quark 
mass as a function of chemical potential (the details are 
shown below in Fig. Q} have been found out, and the 
appearance of state filled up with quarks which is almost 
degenerate with the vacuum state both in quasi- parti- 
cle chemical potential and in ensemble pressure has been 
discovered. 



Here we are going to study the quark ensemble fea- 
tures at finite temperature and fixed baryonic chemical 
potential and to analyse the first order phase transition 
which takes place in such a system of free quasi-particles. 
Analysis is performed within the framework of two ap- 
proaches which are supplementary, in a sense, albeit give 
the identical results. One of those approaches, based on 
the Bogolyubov transformations, is especially informa- 
tive to study the process of filling the Fermi sphere up 
because at this point the density of quark ensemble devel- 
ops a continuous dependence on the Fermi momentum. 
It allows us to reveal an additional structure in the solu- 
tion of the gap equation for dynamical quark mass just in 
the proper interval of parameters characteristic for phase 
transition and to trace its evolution. The result is that 
quark ensemble might be found in two aggregate states, 
gas and liquid, and the chiral condensate is partially re- 
stored in a liquid phase. In order to make these conclu- 
sions easily perceptible we deal with the simplest version 
of the NJL model (with one flavor and one of the standard 
parameter sets) and, actually, are not targeted to adjust 
the result obtained with well-known nuclear liquid-gas 
phase transition. Besides, it seems our approach might 
be treated as a sort of microscopic ground of the conven- 
tional bag model and those states filled up with quarks 
are conceivable as a natural 'construction material' for 
baryons. 

Now as an input to start with we remind the key mo- 
ments of approach developed. The corresponding model 
Hamiltonian includes the interaction term taken in the 
form of a product of two colored currents located in the 
spatial points x and y which are connected by a formfac- 
tor and its density reads as 



H = -q(hV + im)q-qt a lt ,q J dyq't b lv q' {A^l 



>, (1) 
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where q = <?(x), q — q(x), q' = q(y), q' = q(y) are the 



2 



quark and anti-quark operators, 
dp 1 



<?m (x) 



(27r) 3 (2N) I72 [«(PA C )« al (p,v)e ,px + 



+6+(p,s,c)v Qi (p,s,c)e Jpx ] 



(2) 



p\ = — p 2 — to 2 , z-is the colour index, a is the spinor 
index in the coordinate space, a + , a and 6 + , 6 are the 
creation and annihilation operators of quarks and anti- 
quarks, a |0) — 0, b |0) = 0, |0) is the vacuum state of free 
Hamiltonian and m is a current quark mass. The summa- 
tion over indices s and c is meant everywhere, the index 
s describes two spin polarizations of quark and the index 
c plays the similar role for a colour. As usual t a = X a /2 
are the generators of SU(N C ) colour gauge group and to 
is the current quark mass. The Hamiltonian density is 
considered in the Euclidean space and 7 M denote the Her- 
mitian Dirac matrices, fi, v = 1,2,3,4. (A®A'£) stands 
for the formfactor of the following form 



(A a X) = 5 al 



2 G 



N2 _ 1 [H* - y)<W - V ( x - y)] > ( 3 ) 

where the second term is spanned by the relative distance 
vector and the gluon field primed denotes that in the 
spatial point y. The effective Hamiltonian density ([1]) 
results from averaging the ensemble of quarks influenced 
by intensive stochastic gluon field see Ref. [3j. For 
the sake of simplicity we neglect the contribution of the 
second term in ([3]) in what follows. The ground state of 
the system is searched as the Bogolyubov trial function 
composed by the quark- ant i- quark pairs with opposite 
momenta and with vacuum quantum numbers, i.e. 



T|0) 



,(4) 



T = n PiS exp{^[a + (p, s)b + (-p, s) + a(p, s)b(-p, s)]}) 

In this formula and below, in order to simplify the nota- 
tions, we refer to only one complex index which means 
both the spin and colour polarizations. The parameter 
ip(p) which describes the pairing strength is determined 
by the minimum of mean energy 

E=(a\H\a). (5) 

By introducing the 'dressing transformation' we define 
the creation and annihilation operators of quasi-particles 
as A = T a T , B + = T b + T~ l and for fermions 
T -1 = ■ Then the quark field operators are presented 
as 

dp 1 



?(x) 



3(x) 



(2tt) 3 (2b 4 |)i/2 
+B+(p,s) V(p,s) e"**] , 
dp 1 



[A(p,s) U(p,s) 



where U(p, s) = U + (p,s) 74, V(p, s) = V + (p,s) 74 are 
the Dirac conjugated spinors. 

In the paper Ref. [ij the process of filling in the Fermi 
sphere with the quasi-particles of quarks was studied by 
constructing the state of the Sletter determinant type 



\N) 



(7) 



|P|<Pf;S 



which possesses the minimal mean energy over the state 
\N). The polarization indices run over all permissible 
values here and the quark momenta are bounded by the 
limiting Fermi momentum Pp. The momenta and po- 
larizations of states forming the quasi-particle gas are 
marked by the capital letters similar to above formula. 
In all other cases the small letters are used. 

As known the ensemble state at finite temperature T is 
described by the equilibrium statistical operator p. Here 
we use the Bogolyubov-Hartree-Fock approximation in 
which the corresponding statistical operator is presented 
by the following form 



g~ H&pp 



Z = Tr {e-1 3 da PP} , 



(8) 



where some approximating effective Hamiltonian H apv 
is quadratic in the creation and annihilation operators 
of quark and anti-quark quasi-particles A + , A, B + , B 
and is defined in the corresponding Fock space with the 
vacuum state \a) and /3 = T . We don't need to know 
the exact form of this operator henceforth because all 
the quantities of our interest in the Bogolyubov-Hartree- 
Fock approximation are expressed by the corresponding 
averages (a density matrix) 

n(P) = Tr{p A+(P;S)A(P;S)} , 



n(Q)=Tr{p B+(Q; T)B(Q; T)} , 

which are found by solving the following variational prob- 
lem. One needs to determine the statistical operator p in 
such a form in order to have at the fixed mean charge 



Q 4 = Ti{p Q 4 } = V 2N C 



dp 



[n(p) - n(p)} , (9) 



(2tt)3 (2|j P4 |) 1 /2 

+B(p,s)V(p,s) e*P x ] , 
moreover, the transformed spinors U and V are given by 
the following forms 

U (p, s) = cos(» u(p, s) - sm(tp) v(-p, s) , 

V(p, s) = sm(tp) w(-p, s) + cos(ip) v(p, s) . 



[ A + (p, s) U(p, s) e- ipx + where (Q* = -J dx 9*74?) 



, , / [A + (p)A(p) + B{p)B+{p)} , 



(2tt)3 \ Pi \ 



for the diagonal component (of our interest here) and the 
fixed mean entropy 



3 



S = — Tr{pln p} = 
dp 



-V2N C 



(10) 



2tt) 3 



n(p) lnn(p) + (1 — n{p)) ln(l — n(p)) + n(p) ]nn(p) + (1 — n(p)) ln(l — n{p))\ 



(S = — In p) , the minimal value of mean energy of quark 
ensemble 

E = Tr{p H} . 

The definition of mean charge ([9]) is given here up to 
the unessential (infinite) constant coming from permut- 



r 



ing the operators BB + in the charge operator Q4. It may 
not be out of place to remind that the mean charge should 
be treated in some statistical sense because it character- 
izes quark ensemble density and has no colour indices. 

The contribution of free part of Hamiltonian 



Hq = — j dx q(x.) (17V + im) q(x) = 

= J -^^[co^ A+{p-,s)A{p-s) + sme A + {-p■,s)B + {p ] s)- 
+^i\TL9 B(-p;s)A(p;s) -cos6» B(p; s)B + (p; s)] , 



into the mean energy is calculated to be as 
Tr{p H } = J ^3 \P4\ (l-cos0) + 



(2tt)* 
dp 
(2^f 



(11) 



\p4 \ cos 9 [n(j>) + n(p)\ 



where T-Lq = Hq/V2N c is the specific energy and 8 = 2ip. 
It is obvious that a natural regularization has been done 
in the first term of Eq. (fTTj) by subtracting the free 
Hamiltonian contribution H a (without pairing quarks 
and anti-quarks). In our particular case it is quite un- 
affected to have the ensemble energy equal zero at the 
pairing angle equal zero. It explains just an appearance 
of unit in the term containing cos 9. 

The interaction part of Hamiltonian, qt a "f ll qq't ar y l ;q' , 
provides four nontrivial contributions. The term 
Tr{pBB + B' B' + } generates the following components: 
V ai (p, s)t^V^(Q, T)V jk (Q, T)4 7 ^(p, s) 
(adding the similar term but with the changes 
Q, T — ¥ Q' , T" which generates another primed quark cur- 
rent); -2V(Q,T) t a ^V(Q',T')V(Q',T')t b ^V(Q,T). 
Here (as in the following expressions) we omitted all 
colour and spinor indices which are completely iden- 
tical to those of previous matrix element. The term 
Tr{pBAA' + B' + } generates the following nontrivial 
contributions: V(p, s)£YU{(\, t)U(q, t)t b ^V{p, s) 
-V(p, s)t a ^U(P, S)U(P, S)t b ^V(p, s) 



r 



-V(Q, T)t a rU(q, t)U±q, t)t b rV(Q, T) 
+V(Q, T)t a ^U(P, S)U(P, S)t b ^V(Q, T). Aver- 
aging Tr{pAA + A'A' + } gives the following terms: 
U(P,S)t a ^U(p,s)U{p,s)t b ^U(P,S) (adding the 
similar term but with the changes P,S — > P',S'); 
-2U(P, S)t a ^U(P', S')U{P', S')t b ^V(P, S). And 
eventually the nontrivial contribution comes from 
averaging Tr{pA + B + B 1 A'} and it has the form 
V(Q,T)t a YU(P,S)U(P,S)t b ^V(Q,T). All other 
diagonal matrix elements generated by the following 
terms Tr{pAA+B'B'+}, Tr{pBB' + A'+A'} : do not 
contribute at all (their contributions equal to zero). 
Similar to the calculation of matrix elements performed 
at zero temperature in Ref. [J we should carry out 
the integration over the Fermi sphere with the cor- 
responding distribution functions in the quark and 

anti-quark momenta J Pf j0ys — > / J0w[ n {p) + ^(p)] m 
our case of finite temperature. All necessary formulae 
for polarization matrices which contain the traces of 
corresponding spinors could be found in Refs. Q and 
Bearing in mind this fact here we present immediately 
the result for mean energy density per one quark degree 
of freedom as 



»= Wr > £ = E/V > 



where E is a total energy of ensemble, 



4 



w 



G 



dp 

(2^)3 
dp 



(27T) 



\pi\ cos 8[n(p) + n(p)] + 2G 
a (6> - m ) [n(p) + n(p)] 



dp 

(2^)3 ' 
f/q 
(2^)3 



[n(p) + n(p)] 



+ /Aw(l-«**)-G/^Bin(*-* m )y 



</q 

(2^ 



(2tt) 3 
[n(«) + ft(g)] / + 



(12) 



(up to the constant unessential for our consideration 
here). In this formula the following denotes are used 
P — |p|j Q — |q|j 0' = ^(9)1 = ^(P + q) an d the angle 
m (p) is determined by the condition as follows 



sin 9„ 



M 



It is interesting to notice that the existence of such an 
angle stipulates the discontinuity of mean energy func- 
tional mentioned above and found out in [jj. It is quite 
practical to single out the color factor in the four-fermion 

coupling constant as G — Now let us make the fol- 
lowing transformations while integrating in the interac- 
tion terms 

2 j dpf(p) fdq-f dpf(p) J dqf(q) - J dp J dq = 
dpfip) J dq(l - f(q)) - J dp(l - f{p)) J dq, 



and changing the variables p f> q in the last term we 
obtain 



j dp f(p) J dq (1 - f(q)) - J dp J dq (1 - f(q)) 
-J dp (l-/(p)) J dq (!-/(?)) . 



Finally we find out for the mean partial energy it looks 
like 



G 



dp 

dp 

(2^)3 



m\ 



— |p 4 | cos 6 [n{p) + nip) - 1] - 



- Om) [nip) + nip) - 1] 



dq 



(13) 



Mn^-O [n(«)+n(«)-l] / 



We are interested in minimizing the following func- 
tional 

n = E-fiQ 4 -TS, (14) 

where u and T are the Lagrange factors for the chemi- 
cal potential and temperature respectively. The approx- 
imating Hamiltonian which we discussed above H app , is 



constructed just by using the information on E — u Q4 
of presented functional (see, also below). For the specific 
contribution per one quark degree of freedom 

F 

1 2N C ' 1 ' 

we receive 



/ = 



dp 

(2^)3 
dp 



j cos in + n — 1) — u (n — n)] + 



dp 



(2tt) 3 

dp 
(2^r)3 



sin(0-0 m ) (n + n-1) 



dq 

(2^)3 



|P4|" 

9' m ) in' + n'-l)I + 



(15) 



[n Inn + (1 — n) ln(l — n) + n Inn + (1 — n) ln(l — n)] 
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FIG. 1: The residual A for equation (JTTJ) as a function of 
dynamical quark mass M q (MeV) at zero value of temperature 
and the following values of chemical potential n (MeV) — 335 
(the lowest curve), 340, 350, 360 (the top curve). 



Here the primed variables correspond to the momentum 
q. The optimal values of parameters are determined 
by solving the following equation system (df/dO = 0, 
df/dn = 0, df/dn = 0) 



1 



T In ( - - 1 ) = , 



(16) 



- 7 In ( 1-1 

n 



. 



|p 4 | sin0 - M cos (9 - m ) = , 
1^4 cos 9 — /i + M sin (9 — 9 m ) 

\p4 1 cos 9 + /i + M sin (9 — 
where we denoted the induced quark mass as 
M(p) = -2G J^S_(n'+n'-l)sm(9'-9' m )I(p+ q ). (17) 

Turning to the presentation of obtained results in the 
form customary for mean field approximation we intro- 
duce a dynamical quark mass M q parameterized as 



sin(0 - 6 m ) 



|P 4 |-(P 2 + M 9 (P)) 



1/2 



(18) 



and ascertain the interrelation between induced and dy- 
namical quark masses. From the first equation of system 
(fTrJ)) we fix the pairing angle 



sm o = 



p M 



\P4\\Pi\ ' 
and making use the identity 

(|p 4 | 2 - M mf + M 2 p 2 = [p 2 + (M - m) 2 ] \p 4 \ 2 , (19) 
find out that 



cos 9 = ± 



\p 4 \ 2 - m M 

\P4\\Pi\ 



For the clarity we choose the upper sign 'plus'. Then, 
as an analysis of the NJL model teaches, the branch of 



At 300 




Q 4 /3V 

FIG. 2: The chemical potential fi (MeV) as a function of 
charge density Q4 = ( m the units of charge/fm 3 . The fac- 
tor 3 relates the densities of quark and baryon matter. The 
top curve corresponds to the situation of zero temperature. 
The curves following down correspond to the temperature val- 
ues T — 10 MeV, ... , T = 50 MeV with spacing T =10 MeV. 



equation solution for negative dynamical quark mass is 
the most stable one. Let us remember here we are dealing 
with the Euclidean metrics (though it is not a principal 
point) and a quark mass appears in the corresponding 
expressions as an imaginary quantity. Now substituting 
the calculated expressions for the pairing angle into the 
trigonometrical factor expression (obtained in Ref. (3]) 



P 

\P*\ 



— cos ( 



m 

R 



and performing some algebraic transformations of both 
parts of equation we come to the determination 



M q (p) = M(p) 



(20) 



And, in particular, the equation for dynamical quark 
mass (JTTJ) is getting the form characteristic for the mean 
field approximation 



M = -2G 



^3 (n' + n'-l) J(p + q). (21) 



The second and third equations of system (|16[) allow 
us to find for the equilibrium densities of quarks and anti- 
quarks as 



1 



1 



,13 (|p 4 |+m) + 1 



(22) 



and, hence, the thermodynamical properties of our sys- 
tem as well, in particular, the pressure of quark ensemble 



P=- 



dE 
dV ' 



By definition we should calculate this derivative at 
constant mean entropy, dS/dV = 0. This condition al- 
lows us, for example, to calculate the derivative d/i/dV. 



6 



However, this way is not reliable because then the mean 
charge Q4 might change, and it is more practical to intro- 
duce two independent chemical potentials — for quarks 
fi and for anti-quarks fl (following Eq. (I22p with an op- 
posite sign). In fact, it is the only possibility to obey 
both conditions simultaneously. It leads to the following 
definitions of corresponding densities 



1 



1 



e P (l-Ril-z-O + l 



e /3 (|P 4 |+p) + 1 



In fact, this kind of description makes it possible to 
treat even some non-equilibrium states of quark ensem- 
ble (but with losing a covariance similar to the situation 
which takes place in electrodynamics while one deals with 
electron-positron gas). But here we are interested in the 
particular case of p, = fi. The corresponding derivative 
of specific energy ^ might be presented as 



J 



dw 
dV 



dp 



dn dfi 
df~L~dV 



dn dfj, 
djidy 



\P 



4 I cos ( 



2Gsin( 



dq 



sm(6' -6' m )(n' + n' -1)1 




0.00 0.05 0.10 0.15 0.20 

Q 4 /3V 

FIG. 3: The ensemble pressure P (MeV/fm 3 ) as a function of 
charge density Q4 at temperatures T — MeV, ... , T = 50 
MeV with spacing T = 10 MeV. The lowest curve corresponds 
to zero temperature. The dashed curve shows the boundary 
of phase transition liquid-gas, see the text. 



Now expressing the trigonometric factors via dynamical 
quark mass and exploiting Eq. (jTTJ) we have for the en- 
semble pressure 



P 



E 



- V2N r 



dp 

(27T) 5 



dn dfi 
dpdV 



dn dfi 
dfldV 



|P 4 |. (23) 



The requirement of mean charge conservation 



dQ^ Q4 



dV 



V 



-V2N { 



dp 



dn dfi 
df~IdV 



dn dfi 
dfidV 



= 0, (24) 



provides us with an equation which interrelates the 
derivatives dfi/dV, dfl/dV. Apparently, here the regu- 
larized expression for mean charge of quarks and anti- 
quarks is meant Acting in a similar way with the 
requirement of mean entropy conservation, dS /dV = 0, 



we receive another equation as 



dp dn n dfi 



(2tt) 3 dfi 
S 

2N r V 2 



In 



1-ndV 



dp dn 



In 



dfi 



(2tt) 3 dfi 1-ndV 



(25) 



Substituting here T In -r^- 
& 1— n 



1 1 



\Pa\ and T In ■ 



—fi— \ Pi\ after simple calculations (keeping in mind that 
p, = fi) we have with taking into account (l24l) that 



dp 



dn dfi 
dildV 



dndfi' 

d^dv ] |P41 



ST 



{2ir) 3 \dfidV dp.dVJ' 1 2N C V 2 2N C V 2 
Finally it leads for the pressure to the following expres- 



P 



E 
V 



S T 



Qa M 
V 



(26) 



(of course, the thermodynamical potential is ft = —P V). 
At small temperatures the anti-quark contribution is 
negligible, and thermodynamical description can be 
grounded on utilizing one chemical potential fi only. If 
the anti-quark contribution is getting intrinsic the ther- 
modynamical picture becomes complicated due to the 
presence of chemical potential fl with the condition fi = fi 
imposed. In particular, at zero temperature the anti- 
quark contribution is absent and we might receive 



P = -£ 



A 1 Pq 



where fi = y Pp + M 2 (Pp), Pp is the Fermi momentum 

and p q = N/V is the quark ensemble density. 

For clarity, we consider mainly the NJL model [l| in 
this paper, i.e. the correlation function behaves as the 
(5-function in coordinate space. It is well known fact 
in order to have an intelligent result in this model one 
needs to use a regularization cutting the momentum in- 
tegration in Eq. <fTS"|) . We adjust the standard set of 
parameters [f| here with |p| < A, A = 631 MeV, m = 5.5 
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FIG. 4: The fragments of isoterms Fig. [2j|3] see text. Chem- 
ical potential fi (MeV) as a function of pressure P MeV/fm 3 . 
The top curve corresponds to the zero isotherm and follow- 
ing down with spacing 10 MeV till the isotherm 50 MeV (the 
lowest curve). 
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FIG. 5: The dynamical quark mass \M q \ (MeV) as a function 
of chemical potential /i (MeV) at the temperatures T = 
MeV, ... , T = 100 MeV with spacing T = 10 MeV. The most 
right curve corresponds to zero temperature. 



MeV and GA 2 /(2ir 2 ) = 1.3. This set of parameters at 
n = 0, n = 0, T = gives for the dynamical quark mass 
M q = 335 MeV. In particular, it may be shown the fol- 
lowing representation of ensemble energy is valid at the 
extremals of functional (fT5")) 

E = E vac + 2N C V { |P 4 | ( n + n) , 

A M2 (27) 

E vac = 2N c Vj j^(\pA - \PA) + 2N - V ^q > 

It is easy to understand this expression with the vacuum 
contribution subtracted looks like the energy of a gas of 
relativistic particles and anti-particles with the mass M q , 
and coincides identically with that calculated in the mean 
field approximation. 

Let us summarize the results of this exercise. So, we 
determine the density of quark n and anti-quark n quasi- 
particles at given parameters /i and T from the second 
and third equations of system (1161) . From the first equa- 
tion we receive the angle of quark and anti-quark pair- 
ing 6 as a function of dynamical quark mass M q which is 
handled as a parameter. Then at small temperatures, be- 
low 50 MeV, and chemical potentials of dynamical quark 
mass order, /i ~ M q , there are several branches of solu- 
tions of the gap equation. Fig. Q] displays the difference 
of right and left sides of Eq. (fTT)) which is denoted by A 
at zero temperature and several values of chemical po- 
tential fj, (MeV) = 335 (the lowest curve), 340, 350, 360 
(the top curve) as a function of parameter M q . The zeros 
of function A(M q ) correspond to the equilibrium values 
of dynamical quark mass. 

The evolution of chemical potential as a function of 
charge density Q4 = & (in the units of charge//m 3 ) 
with the temperature increasing is depicted in Fig. [2] 
(factor 3 relates the quark and baryon matter densities) . 
The top curve corresponds to the zero temperature. The 
other curves following down have been calculated for the 



temperatures T = 10 MeV, ... , T = 50 MeV with spac- 
ing T = 10 MeV. As it was found in Ref. [4| the chemi- 
cal potential at zero temperature is increasing first with 
the charge density increasing, reaches its maximal value, 
then decreases and at the densities of order of normal 
nuclear matter density^, p q ~ 0.16//m 3 , becomes al- 
most equal its vacuum value. Such a behaviour of chem- 
ical potential results from the fast decrease of dynamical 
quark mass with the Fermi momentum increasing. It is 
clear from this Fig. the charge density is still multivalued 
function of chemical potential at the temperature slightly 
below 50 MeV. The Fig. [3] shows the ensemble pressure 
P (MeV/fm 3 ) as the function of charge density Q4 at 
several temperatures. The lowest curve corresponds to 
the zero temperature. The other curve following up cor- 
respond to the temperatures T = 10 MeV, ... , T = 50 
MeV (the top curve) with spacing T = 10 MeV. It is 
interesting to remember now that in Ref. Q the vac- 
uum pressure estimate for the NJL model was received 
as 40—50 MeV/fm 3 which is entirely compatible with the 
results of conventional bag model. Besides, some hints at 
instability presence (rooted in the anomalous behavior of 
pressure dP/dn < 0) in an interval of the Fermi momenta 
has been found out. 

Fig. 0]shows the fragments of isotherms of Fig. [21[3]but 
in the different coordinates (chemical potential — ensem- 
ble pressure) . The top curve is calculated at the zero tem- 
perature, the other isotherms following down correspond 
to the temperatures increasing with spacing 10 MeV. The 
lowest curve is calculated at the temperature 50 MeV. 
This Fig. obviously demonstrate a presence of the states 
on isotherm which are thermodynamically equilibrated 
and have equal pressure and chemical potential (see the 
characteristic Van der Waals triangle with the crossing 
curves). The calculated equilibrium points are shown 
in Fig. [3] by the dashed curve. The intersection points 
of dashed curve with an isotherm are fixing the bound- 
ary of gas — liquid phase transition. The corresponding 
straight line P = const which obeys the Maxwell rule 
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FIG. 6: The dynamical quark mass \M q \ (MeV) as a function 
of temperature at small value of charge density Q4 . 

separates the noii-equlibrium and unstable fragments of 
isotherm and describes a mixed phase. The correspond- 
ing critical temperature for the parameter we are using 
in this paper turns out to be T c ~ 45.7 MeV with the 
critical charge density as Q4 ~ 0.12 charge/fm 3 . Usu- 
ally the thermodynamic description is grounded on the 
mean energy functional which is the homogeneous func- 
tion of particle number like E = N f(S/N, V/N) (with- 
out vacuum contribution). It is clear such a description 
requires the corresponding subtractions to be introduced, 
however, this operation does not change the final results 
considerably. It was argued in Refs. [4] that the states 
filled up with quarks and separated from the instability 
region look like 'natural construction material' to form 
the baryons and to understand the existing fact of equi- 
librium between vacuum and octet of stable (in strong 
interaction) baryons [loj. 

The dynamical quark mass \M q \ (MeV) as a function 
of chemical potential \x (MeV) is presented for the tem- 
peratures T = MeV, ... , T = 100 MeV with spacing 
T = 10 MeV in Fig. [5j The most right-hand curve corre- 
sponds to the zero temperature. At small temperatures, 
below 50 MeV, the dynamical quark mass is the multi- 
valued function of chemical potential. The Fig. [5] shows 
the dynamical quark mass as a function of temperature 
at small values of charge density Q4 ~ 0. Such a be- 
haviour allows us to conclude that the quasi-particle size 
is getting larger with temperature increasing. It becomes 
clear if we remember that the momentum corresponding 
the maximal attraction between quark and anti-quark pg 
(according to Ref. Q) is denned by dsm9/dp — 0. In 
particular, this parameter in the NJL model equals to 

pg = (\M q \ to) 1 / 2 . (28) 

but its inverse magnitude defines the characteristic (ef- 



fective) size of quasi- particle rg = p g . 

If one is going to define the quark chemical potential as 
an energy necessary to add (to remove) one quasi-particle 
(as it was shown in Q at zero temperature), \i — dE/dN, 
then in vacuum (i.e. at quark density p q going to zero) 
quark chemical potential magnitude coincides with the 
quark dynamical mass. It results in the phase diagram 
displayed at this value of chemical potential although, 
in principle, this value could be smaller than dynamical 
quark mass as it has been considered in the pioneering 
paper [6|. If one takes, for example, chemical potential 
value equal to zero it leads to the conventional picture 
but, obviously, such a configuration does not correspond 
to the real process of filling up the Fermi sphere with 
quarks. 

Apparently, our study of the quark ensemble thermo- 
dynamics produces quite reasonable arguments to pro- 
pound the hypothesis that the phase transition of chiral 
symmetry (partial) restoration has already realized as the 
mixed phase of physical vacuum and baryonic matter [lll|. 
However, it is clear our quantative estimates should not 
be taken as ones to be compared with, for example, the 
critical temperature of nuclear matter which has been 
experimentally measured and is equal to 15 - 20 Mev. 
Besides, the gas component (at T = 0) has non-zero den- 
sity (as 0.01 of the normal nuclear density) but in reality 
this branch should correspond to the physical vacuum, 
i.e. zero baryonic density [l^. In principle, an idea of 
global equilibrium of gas and liquid phases makes it pos- 
sible to formulate the adequate boundary conditions at 
describing the transitional layer arising between the vac- 
uum and filled state and to calculate the surface ten- 
sion effects. We are planning to consider these aspects of 
phase transition in a separate paper. 

As a conclusion we would like to emphasize that in the 
present paper we demonstrated how a phase transition 
of liquid-gas kind (with the resonable values of parame- 
ters) emerges in the NJL- type models in which the quarks 
are considered as the quasi-particles of the model Hamil- 
tonian with four-fermion interaction. The constructed 
quark ensemble displays some interesting features for the 
nuclear ground state (for example, an existence of the 
state degenerate with the vacuum one) but needs further 
study of its role in the context of existing research {8j 
activity to explore the complicated (or, may be, more re- 
alistic) versions of the NJL model and knowledge of the 
QCD thermodynamics as gained in the lattice simula- 
tions. 

The authors are deeply indebted K. A. Bugaev, R. N. 
Faustov, S. B. Gerasimov, K. G. Klimenko, E. A. Ku- 
raev, A.V. Leonidov, V. A. Petrov and A. M. Snigircv 
for numerous fruitful discussions. 
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